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Abstract 



Model calculations are performed of extensive air shower (EAS) component en- 
ergies using a variety of hadronic interaction parameters. A conversion factor from 
electromagnetic component energy to the energy of ultra-high energy cosmic rays 
I> ' (UHECRs) and its model and primary mass dependence is studied. It is shown 

that model dependence of the factor minimizes under the necessary condition of the 
■^-j- \ same maximum position and muon content of simulated showers. 

CO 

: 1 Introduction 

o 

UHECR particles hitting Earth atmosphere produce a cascade of secondary particles, 
the small part of which is detected on the ground with EAS array. The energy of the 
primary particle, Eq, is distributed among the shower components. The most of the energy 
deposit is due to ionization and excitation of the air molecules caused by electromagnetic 
component. 

The primary energy estimation algorithms are based mainly on the measurement of the 
shower parameters related to electromagnetic component energy, E em . Such an approach 
is realized in High Resolution Fly's Eye (HiRes) pQ, Pierre Auger Observatory (PAO) [2] 
and the Yakutsk array [31 H] experiments. Future applications are planned in the Telescope 
Array and satellite projects. 

In all these experiments a fraction of the primary energy cannot be measured because it 
is carried away by hadrons, muons etc., unobservable with array detectors. This 'missing 
energy' and conversion factor E em /E can be calculated modeling a cascade in atmosphere, 
as was done in [UE1E1E]- However, only QGSJET and SIBYLL models have been applied 
in these simulations, setting aside the model dependence of the conversion factor. 
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In the present work the energies of air shower components are calculated using ex- 
tremely different hadronic interaction models, but having the same maximum position 
and muon content of the resultant showers. It is shown that the model dependence of the 
ratio E em /E minimizes in this case. The variability of the ratio is estimated due to mod- 
els and the primary nucleus mass. The experimental uncertainties in the electromagnetic 
component energy measurements are discussed. In particular, the energy spectra mea- 
sured with the Yakutsk and HiRes arrays are shown to be coincident within experimental 
errors. 

2 Basic experimental data 

Relativistic electrons of the shower induce Vavilov-Cherenkov radiation (Cherenkov light) 
in the atmosphere. The total flux of light, Qtot, can be used as an estimator of the 
electromagnetic component energy in EAS. This method realized in the Yakutsk array 
experiment is based on the relation between Q tot and E em largely independent of the 
model [3l [8] . The energy fraction dissipated in the ground is estimated using the tail of 
cascade curve measurement in inclined showers. 

Another method, namely fluorescence technique, is applied in HiRes and PAO/FD 
experiments where the fluorescence light emitted by wounded nitrogen molecules along 
the trajectory of the shower is collected by mirrors and received by photomultiplier pixels. 
The light intensity is proportional to the number of shower electrons, so the cascade curve, 
N e (t), is scanned by fired pixels and ionization in the atmosphere (ionization integral) can 
be estimated: 



where t is slant depth along the shower axis in units of radiation length (36.7 g/cm 2 [9]); 
to is the thickness of atmosphere; energies are in units of critical energy (eo = 86 MeV). 
As it was shown in cascade theory [10], the electromagnetic component energy is equal to 
ionization integral E em = E?°; that is why the energy E em is actually evaluated in these 
experiments. 

Measurements of the shower maximum position in the atmosphere as a function of the 
primary particle energy, t max (E ), and the number of muons at observation level (usually 
as a ratio to the number of electrons, N^/Ne) are the basic parameters to be guided by 
performing the cascade simulations in the UHE region. The only data available in the 
region are those given by the HiRes, AGASA and Yakutsk experiments [Tj, [31 [EE]. In the 
HiRes case the position of the shower maximum is measured by stereo system of two Eyes. 
On the contrary, the Yakutsk array detectors are measuring the Cherenkov light, and the 
lateral distribution form of the light spot on the ground is used to derive t max in this case 



The number of electrons and muons at the ground level is measured by scintillation 
detectors/proportional counters (shielded with Fe/concrete and ground to detect muons) 
of the Yakutsk and AGASA arrays (Fig. ED. 




(Fig. HJ . 
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Figure 1: Shower maximum depth in the atmosphere as a function of EAS primary particle energy. 
HiRes pQ (triangles) and HiRes/MIA [T^] (rhombuses) results are shown in comparison with the Yakutsk 
array data (circles). 

Experimental errors in t max measurement are 0.82 to 1.38 in the energy range Eq€q G 
(10 17 , 10 19 ) eV; N^/N e ratio differs up to ~ 2 times between Yakutsk and AGASA data 
while the difference in iV M itself is around 25% [TT]. The main source of the discrepancy is 
the total number of charged particles on the ground which is estimated using the particle 
density measured at the shower periphery in both experiments. 

3 Modeling the cascade in atmosphere 

As was shown in previous calculation of the energy balance of EAS components [I], 
the dependence of primary energy deposited in the atmosphere on hadronic interaction 
models (i.e. multiplicity of secondaries, cross sections, etc.) can be parameterized via the 
longitudinal shower development characteristics t max and N^/N e . This is a consequence of 
the close connection between E em and the ionization integral. In this case (as far as other 
shower parameters are not concerned), we can go far beyond the bounds of conventional 
models. 

Namely, we can use in very high energy region a variety of arbitrary extrapolations 
of accelerator data in order to apply hadronic interaction models with only restriction of 
the same resultant shower maximum and the muon content. It will allow to estimate the 
real limits of the ratio E em /E variations using distinctly different models in addition to 
'standard' models fitted to describe the same experimental data. 

In this work three hadronic interaction models which are characterized by extremely 
different forms of the rapidity, y, distribution of the 'sea' secondaries in pionization region 
are used: i) Gauss model with normal rapidity distribution; ii) Delta model with equal 
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Figure 2: The number of muons (E^ > 1 GcV) vs the number of electrons at the depth to = 27.8. 
Current and previous [3] Yakutsk array data are shown by squares and triangles. AGASA data approxi- 
mation [llj is given by the dash-and-dot line. 

rapidities of n s secondaries and iii) Flat distribution model with f(y) = const. All other 
parameters of the models are flexible in order to get t max and N^/N e of the shower 
fitting experimental data within errors. In the Flat distribution model, for instance, 
the pionization region width and the fraction of secondary pions, kaons and nucleons 
in multiple production processes have been adjusted in addition to cross sections and 
fragmentation coefficients to get the output EAS observables desired. 

Numerical solution of hadron transport equations is obtained using the resolvent of 
Volterra equation on the rectangular lattice {ij, y^} [13] . To calculate the number of 
electrons in the shower at a depth t, Greisen's formula is employed. Electromagnetic 
component energy is assumed to be equal to the energy of neutral pions, the decay of 
which into gamma-quanta initiates the component. The muon and neutrino energy is 
calculated via the energy fraction taken from the decay of charged pions and kaons. The 
hadronic component energy is summed up at each layer t^. 

It is well-known that t max is strongly fluctuating parameter of the shower. Monte Carlo 
simulations using CORSIKA code give RMS deviation 1.9 and 1.1 for the proton and Fe 
initiated showers of the same Eq, correspondingly [H]. On the contrary, the primary 
energy fraction deposited in atmosphere is fluctuating much lesser because it is related to 
the integral along the shower trajectory. 

In order to estimate the influence of cascade fluctuations on the ionization integral, 
a leading fragment approximation in the simplest case of the primary nucleon is used 
- only fluctuations in inelasticity and interaction points of the EAS primary particle 
are considered, along suggestions given in [TS]. Bearing in mind the isotopic invariance 
of multiple production processes and multiplicity distribution, we derive the ionization 
integral event-by-event distribution, f(E em ) with fixed E , as a {n = xq/X ~ 10 - fold) 
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Figure 3: Event-by-event fluctuations in a leading fragment approximation. 

convolution of inelasticity fluctuations. Resultant distribution is a Gaussian with average 
energy fraction 

E^> 1 - K n 



E l-K 

where K is an elasticity in nucleon interactions; K n o is a projectile energy fraction carried 
out by neutral secondary pions. The RMS deviation is ~ 0.1% as illustrated in Fig. [31 We 
neglect fluctuations in the rest of Eq distributed among the shower components because 
of averaging over the huge number of shower particles. 



4 Results 

We are interested in the ionization integral and the energy deposited to muons+neutrinos 
in a shower. These values calculated for our models (Eoe c = 10 18 eV) are presented in 
Fig. HI The hadronic component energy is less than O.OlE'o at sea level. 

Ionization in the atmosphere, E° } is closely related to the total air Cherenkov light 
flux on the ground [31 H] and is measured with Cherenkov light detectors of the Yakutsk 
array. A point in Fig. H] is given with experimental error in the aggregate of calibration, 
atmospheric extinction of light uncertainty, etc. 

The energy of muons and neutrinos is estimated as a sum of measured on the 
ground and model calculation results giving the energy of neutrinos, ionization and decay 
of muons [T7] . Absolute experimental E^+ v uncertainty (~ 2%) is based on the difference 
between Yakutsk and AGASA measurements. 

Another presentation of the results is given in Fig. [5j Electromagnetic component 
energy is E em = E\° + Ei, where Ei is the energy of electrons and photons dissipated in 
the ground. To estimate it we need the total number of electrons at t and its attenuation 
length: E\ ~ \N e N e (to). The result based on the Yakutsk array data is shown in Fig. [5] 
by circles. 
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Figure 4: Energy fractions of EAS components as a function of the depth in the atmosphere. Model 
calculation results are given for hadronic (h), electromagnetic (e+7) components and muons and neutrinos 
{[i + v). The solid line refers to Gauss model, dotted line to Delta model and dash-and-dot line to the Flat 
rapidity distribution model. Experimental data of the Yakutsk array are shown at t — 27.8: ionization 
in the shower (circle) and the energy of muons+neutrinos (square) [IS] . 

The Yakutsk array data are used to estimate the energy of muonic component. Energy 
remainder deposited to undetectable hadrons and neutrinos is estimated using the cascade 
modeling. In the HiRes case, the energy fraction E em /E is exclusively the product of 
model simulations because the only component measured is the fluorescence light. Since 
CORSIKA/QGSJET code is used to estimate the conversion factor to primary energy, 
triangles represent QGSJET model results. EAS simulation with SIBYLL model gives 
the ratio with difference less than 1.6% [6]. 

The main feature of results given in Figs. B1I51 is that different models lead to ap- 
proximate fractions of the primary energy assigned to electromagnetic component, if 
t max and N^/N e are coincident. The model uncertainty is less than 3% in the interval 
E^eo G (10 17 , 10 19 ) eV where the experimental data are available; the energy dependence 
of E em / E is limited to ~ 2± 1% per decade, while a spread at 10 18 eV is within 87 ± 1%. 

Additionally, the ionization integral and E^+v are varying with the primary particle 
mass, in spite of the same t max for different nuclei. We estimated the variance to be 
less than 6% if the primary mass is 1 < A < 56. The sum of two uncertainties (from 
the model and primary mass) is 5E em /Eo < 7%. We have not considered here showers 
initiated by other primaries, such as photons, neutrinos or mini black holes, because their 
muon content is a subject of the alternate investigation. 

Due to the negligible fraction of hadronic component energy at sea level, the remainder 
of E is transferred to muons and neutrinos. As a consequence, E^ +u /E is nearly model 
independent, too, but the muon component energy measurable on the ground is more 
variable in different models. 
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Figure 5: Energy fraction of electromagnetic component vs the primary particle energy. The hadronic 
interaction model results are shown by the same signs as in Fig. A HiRes estimation is shown by 
triangles pQ. Simulation uncertainty 6% due to unknown primary nucleus mass is assigned to points. 
Two circles at EqEq = 10 18 and 10 19 eV illustrate E em /Eo evaluation based on the Yakutsk array 
data [4]. Experimental errors are shown by vertical bars. 



5 Experimental uncertainties 

The experimental errors in the shower maximum detection and the number of muons 
at to lead to the uncertainty of calculated ratio E em / E (we will pass over SE^/Eq 
considering it congruous). We have modeled this varying the multiplicity of secondaries 
and cross sections which result in the changes of t max and N^E > lGeV) comparable to 
experimental uncertainties. Resultant SE em /E turned out to be below 5% in the case 
of depth variation 5t max < 1.36, and 3% for SN^/N^ < 0.25. So we have assumed the 
aggregate experimental uncertainty due to t max and measurement errors to be below 
6%. 

The greatest uncertainty source is the ionization integral measurement itself. For the 
Yakutsk array data it comprises of errors due to Cherenkov light measurement and the 
number of charged particles on the ground |H[8j: uncertainty in atmospheric transparency 
(15%); detector calibration (21%) and total light flux measurement (15%) errors; an 
uncertainty in the number of electrons reaches 60% because the only measurable parameter 
is the particle density beyond hundred meters from the shower core. Resultant E em / E 
estimation uncertainty (~ 30%) is the sum of all these errors weighed with the shower 
component energies. 

The HiRes group claims a systematic uncertainty ~ 20% [1] aggregated of errors in 
the absolute calibration of the photo-tubes (10%), the yield of the fluorescence process 
(10%), modeling of the atmosphere, and so on. 

However, a comparison of model predictions with the HiRes/MIA measurement of the 
average cascade curve in the interval (10 17 , 10 18 ) eV (Fig. [6]) reveals a discrepancy between 
calculated and measured curves. Substantially, a contradiction in the cascade curve can be 
resolved using the attenuation length of electrons measured at sea level. Equal intensity 
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Figure 6: The average longitudinal profile of EAS (circles,i?oeo S (fO 17 , 10 18 ) eV) measured in the 
HiRes/MIA hybrid experiment [T2] and shower cascade curves calculated in models. Signs for models 
are the same as in previous Figures. The average cascade curve of electrons (triangles) for proton- 
initiated vertical showers of energy EqCq — 10 175 eV simulated using CORSIKA/QGSJET [16] is given 
for comparison. Curves are adjusted to t max = 19.1. 

cuts method gives the value Ajv e = 5.6 ± 0.9 in the interval t G (27.25,29.97) for the 
Yakutsk array data [3J, while the HiRes curve has Ajv e = 15.9 ± 1.1. Models predict the 
length Xn £ £ (4.7,5.4) in agreement with the Yakutsk array data. 

There may be an additional uncertainty up to 40% in E em assignment due to inad- 
equate fitting of the longitudinal shower profile detected with HiRes. This may be the 
result of increased measurement error of the faint fluorescence light far from the shower 
maximum, and/or the direct and scattered Cherenkov light contamination to the signal. 

There is an alternative E em estimation method proposed in [7]. It is based on the 
relation between the number of electrons at the shower maximum, N max , and energy of 
gamma-quanta initiating electromagnetic sub-cascade: N max = 0.3E 1 / JhxE 1 [10J. An 
advantage of the method is in avoidance of the cascade curve approximation; disadvantage 
consists in the asymmetric distribution of N max . On the whole, it seems to be the more 
reliable approach to the estimation of electromagnetic component energy in comparison 
with integration of the longitudinal shower profile, not only at energies above 10 19 eV, as 
was concluded in [7]. 

Presumably, the overall experimental uncertainty in ionization integral estimation is 
~ 30% for the Yakutsk array data, and < 45% for the HiRes measurements. Consequently, 
notwithstanding nearly the same conversion factor from E em to Eq, the primary energy 
estimate may be distinctly different in these experiments due to systematic errors in E em 
evaluation. 

To provide an illustration, the differential energy spectra measured with the Yakutsk 
array [4] and HiRes pQ are shown in Fig. [7j The primary energy is decreased by 20% for 
the former data, and increased by 30% for the latter. Intensities are multiplied by i?o e o 
in order to demonstrate the spectrum features - 'ankle' around 10 19 eV and steepening at 
E e ~ 10 20 eV (Greisen-Zatsepin-Kuzmin effect). The shape and intensity of the spectra 
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Figure 7: The energy spectrum of cosmic rays. Monocular data from two HiRes Eyes are given: HiRes- 
I (filled triangles) and HiRes-II (open triangles). The Yakutsk array data are selected with trigger-500 
(open circles), and trigger- 1000 (filled circles). Energy correction factors are given in the text. 



derived from the data of the Yakutsk array and HiRes are consistent indicating a need for 
systematic correction (within experimental errors) to the ionization integral measured. 



6 Conclusion 

Extremely different hadronic interaction models lead to the resultant (electromagnetic 
component energy/EAS primary particle energy) ratio consistent within 7% if the model 
and primary mass guarantee the same maximum position in atmosphere and muon content 
of the shower: 

E Ea 

= (0.87 ± 0.01) + (0.02 ± 0.01) lg 
£■0 10 

10 17 < E e < 10 19 eV. 

The primary energy estimation algorithms based on ionization integral measurement 
rely on this ratio and the resultant uncertainty originates from experimental errors pre- 
dominantly; model dependent one is minor in the presence of contemporary measurements 
of the longitudinal shower development parameters. 
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